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Abstract. We have performed a heavy ion radiative capture reaction between two Ught heavy ions, '^C and 
'*0, leading to -'^Si. The present experiment has been performed below Coulomb barrier energies in order to 
reduce the phase space and to try to shed light on structural effects. Obtained y-spectra display a previously 
unobserved strong feeding of intermediate states around 1 1 MeV at these energies. This new decay branch is 
not fully reproduced by statistical nor semi-statistical decay scenarii and may imply structural effects. Radiative 
capture cross-sections are extracted from the data. 



1 Introduction 

The isotope ^^Si can be considered as a key nucleus to un- 
derstand the coexistence between mean field effects and 
cluster structures. Indeed different structures coexist in its 
excited states even at low excitation energies. Recent anti- 
symmetrized molecular dynamics (AMD) calculations for 
^^Si LU have shown that cluster structures such as '^C- 
'^O or ^''Mg-Q' may have large contributions for normal- 
deformed and superdeformed states, respectively. By us- 
ing the heavy-ion radiative capture (HIRC) mechanism in 
which the compound nucleus decays solely by y-ray emis- 
sion we have investigated the '^C-"'0 cluster structure in 
^^Si. Moreover this mechanism will directly populate ^^Si 
at energies where the calculated nuclear matter densities 
have similar asymmetries as the '^C-h'^O reaction [2]. 

The excitation function of '^C-h'^O exhibits naiTow res- 
onances |3| around the Coulomb BaiTier (Vb ~ 7.8 MeV). 
These resonances are correlated in all reaction channels in- 
cluding the HIRC. Baye and Descouvemont |4 1 have shown, 
using GCM calculations, that for the similar system '^C-H 
'^C, several observed resonances can be interpreted in terms 
of molecular resonances. Although HIRC is a rare process 
compared to the dominant fusion-evaporation channels (n, 
p, a), y-decay will favourably populate states with large 
structural overlaps with the entrance-channel thus allow- 
ing us to try to identify the possible cluster states. 



2 Experiment 

Our experimental program is focused on the study of HIRC 
in the '^C-h'^c and '^C-h'^'O ^M systems. The present 
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'^C('^0,y) reaction has been performed at the TRIUMF 
facility (Vancouver, Canada). A '^O ISAC beam has been 
used at two resonant energies down to 15% below Vb (Eiab 
= 1.07 and 0.96 AMeV) on high purity (99.9%) '^c thin 
(50 fj.g.cm^^) targets. The DRAGON recoil separator has 
been used at 0° to select the ^^'Si recoils. We have used 
the BGO array to record the y-rays in coincidence with the 
heavy recoils identified in a DSSSD placed at the DRAGON 
focal plane. 



3 Results and discussion 

3.1 y-ray decay 

Fig. [T] displays the highest energy y-ray in each y-event 
recorded in the BGO array in coincidence with the ^^Si 
for the two studied energies. Both spectra show the same 
global structure. Concerning the region above 17 MeV, the 
direct transition to the ground state (g.s.) contributes to the 
bump above 20 MeV but due to the response function of 
the spectrometer and the deviation out of 0° of the recoils 
induced by high energy y-ray emission very few recoils en- 
ter DRAGON. The feeding of the 2| ( 1 .778 MeV) from the 
resonances is also observed (Ey^ ~ 21 MeV) at both ener- 
gies. For the lowest energy we also have a y-line around 18 
MeV corresponding to the feeding of the 4| (4.617 MeV). 
Concerning the low energy part of the spectra, the de- 
cay of the two first excited states of the ^**Si can be identi- 
fied {Ey„ = 1.78 and 2.84 MeV). The decay of the prolate 
head band O3 (6.690 MeV) to the 2| contributes to the line 
around 5 MeV but this line does not correspond to a sin- 
gle state decay. We observe also the decay of the 3 j^ to the 
g.s. which contributes to the y-line at 6.8 MeV, a direct 
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Fig. 1. Highest energy y-ray spectra in coincidence witli -'*Si for 
tlie '-C+'*0 reaction at £„„ = 6.6 MeV (daslied line) and Ec_,„ = 
7.2 MeV (full line). Spectra are normalized to the same integral. 



decay from the resonance is also observed for the two en- 
ergies. It has been shown at higher energies [51 that this 
first state of negative parity in ^'^Si is crucial to understand 
the y-spectra. 

Due to the low resolution of the BGO y-aiTay and the 
large number of states in ^'^Si around 10 MeV, the large 
bump between 10 and 15 MeV with correspond to the feed- 
ing of intermediate states around 8-13 MeV cannot be dis- 
cussed in terms of the feeding of particular states. More- 
over in this region the transitions from the resonances and 
the subsequent decays to the g.s. are quite close in en- 
ergy. This previously unobserved y-flux at these reaction 
energies has to be linked to what has been observed at 
higher energies (Ec,i„ - 8.5, 8.8 and 9.0 MeV) |5J, where 
the bump was around 14 MeV. Taking into account that 
in our previous experimental campaign the ^^Si was pop- 
ulated at higher excitation energies {E* ~ 25.3, 25.6 and 

25.8 MeV), than in the present experiment {E* ~ 23.3 and 

23.9 MeV), the observed differences on the bump centroids 
are only resulting from the differences in entrance-channel 
energies. 
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Fig. 2. Spin distributions normalized to fusion cross-sections ob- 
tained by CCFULL |8| calculations for the two studied energies: 
red triangles correspond to E^j,, = 6.6 MeV and blue squares to 
£,,,, = 7.2 MeV. 
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Fig. 3. Highest energy y-ray spectra in coincidence with -'*Si for 
£■,..„,= 6.6 MeV (full line) £„„.= 7.2 MeV (full line) compared 
to fully statistical decay (red dashed line) using the spin distri- 
butions given in Fig.|2] Simulation spectra are normalized to the 
data integrals. 



In order to fully understand the decay pattern of the 
resonances and see if we have a structural effect which 
will favourably feed particular states, we compared our -y- 
spectra to GEANT simulations with different conditions 
in the entrance-channel. All simulations include 68 known 
bound or quasi bound states of ^**Si between and 13 MeV 
and their y-decays. Above the particle threshold (~ 10 MeV 
for the lowest one, in this case Q'-i-^'*Mg) we have selected 
states with a large Fy compared to particle emission. To es- 
timate the branching ratios of the entrance-channel to each 
of these states we have used Weisskopf estimates and the 
reported average strengths of electric and magnetic transi- 
tions in ^^Si |7|. In this self-conjugate nucleus (T. - 0), 
some particular selection rules on isospin apply for L=l, 
AT-Q transitions: El are forbidden and Ml strengths are 
reduced by a factor of 100. As we use an entrance-channel 
with r = in all simulations, these rules will strongly in- 
fluence the obtained y-spectra. 

We have first simulated a fully statistical scenario, which 
consists of a spin distribution in the entrance-channel. Spin 
distributions for the two studied energies are given in Fig.|2] 
These distributions are the results of coupled-channel cal- 
culations, performed with the CCFULL code |I8], and are 
obtained by adjusting the diffuseness parameter (a =0.57 
and 0.55 for /if ,„ =6.6 and 7.2 MeV) in order to reproduce 
the fusion cross-sections from [ 1 1 1 and given in Fig.|2] The 
two curves are normalized to the calculated fusion cross- 
sections. Below Vb, phase space is reduced and spin dis- 
tributions are quite narrow and centered at low spin (2^) 
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compared to what was obtained at higher energies for the 
i2c+"'0 reaction |5|. 

Simulated spectra (dashed lines) are given in Fig. [3] 
along with experimental data (full lines). For the two ener- 
gies we globally reproduce the y-spectrum shape. Indeed a 
large bump between 11 and 15 MeV dominates the spec- 
trum. However looking into the details we see that this sta- 
tistical scenario is unable to reproduce the intensities of the 
different transitions between the low-lying states. Another 
discrepancy concerns the high energy part of the simula- 
tions where the direct feeding of the two first excited states 
is overestimated. Concerning the intermediate region of the 
spectra the bump is better reproduced for the highest en- 
ergy. For the lowest energy, the fully statistical model has 
too large branching ratios to states lying below 10 MeV. 
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Fig. 4. Highest energy y-ray spectrum for E^,,, = 6.6 MeV (full 
line) compared to decay scenarii with two different spins in the 
entrance-channel: 1" (red dashed line) and 2* (red dotted line). 
Simulations spectra are normalized to the data integral. 



As the studied energies correspond to resonant ener- 
gies, we used another scenario which consists of a unique 
spin in the entrance-channel. As the spin distributions given 
in Fig. |2] show that entrance-channel spins greater than 5H 
have very small contributions to the fusion cross-sections, 
we limit the discussion here to entrance spins less or equal 
fo 4S. Fig.|4]displays the simulation results for resonances 
J" = 1^ and 2^ and their comparison with the lowest en- 
ergy data. As for the first scenario, we see that we overall 
reproduce the shape of the spectrum with the negative par- 
ity entrance spins, even if the bump seems to be shifted. 
Again this is due to too large branching ratios to states ly- 
ing below the fed states in the experiment and has to be 
interpreted in terms of structural effects. In the case of a 
positive parity resonance the large y-flux going directly to 
the g.s. and the 2| depletes the region around Ey„ = 12 
MeV and gives rise fo a completely different y-spectrum. 

In order to have a objective criterion to compare simu- 
lated and experimental spectra we used the so-called ^^ test 
of homogeneity. Evolution of the x^ versus the resonance 
spin for the two studied energies is given in Fig.|5] Looking 
at the two curves on Fig. |5] we observe a kind of oscilla- 
tion between even and odd spins. This can be understood 
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Fig. 5. ^- variation as a function of the spin of the entrance- 
channel, red triangles correspond to Ec.m. = 6.6 MeV and blue 
squares to £■„„. = 7.2 MeV. 



by the fact that for the even spins, corresponding to posi- 
tive parity in the entrance-channel, decays to the low-lying 
states of positive parity are mainly achieved by E2 transi- 
tions which are not slowed down by any selection rules in 
this nucleus. 




Fig. 6. Highest energy y-ray spectrum for E^,,, = 7.2 MeV (full 
line) compared to decay scenarii with two different spins in the 
entrance-channel: 3" (red dashed line) and 2^ (red dotted line). 
Simulations spectra are normalized to the data integral. 



For the highest energy for which results are given in 
Fig.|6]the same discussion can be held. The difference is 
that a y"^ = 3" resonance in the entrance-channel repro- 
duces the data better than a J" = 1". But as for the lowest 
energy, data are better reproduced with entrance-channel 
odd spins. Branching ratios to negative parity states from 
an negative entrance spin are expected to be favored against 
branching ratios to positive parity states in resonance y- 
decay. This is partially reproduced by the semi-statistical 
scenario. We already stressed the importance of particu- 
lar negative states such as the 3j^ which is the first nega- 
tive parity state, at higher energies. As the phase space is 
reduced due to the reduction of the bombarding energy. 
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this particular decay may be more important in the two 
explored energies than in our previous experimental cam- 
paign. 



3.2 Radiative capture cross-section 

The previously unobserved large feeding of intermediate 
states around 1 1 MeV allows us to reevaluate the radiative 
capture (RC) cross-section measured in Ref. flQ\. 
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in which the cross-section, cri^c, is given in barn. A^^, A^, 
and A^, correspond, respectively, to the number of ^^Si re- 
coils observed at DRAGON focal plane (±10%), the num- 
ber of '^C atoms/cm^ (±10%) in the target and the number 
of incident "'O per second (±10%). Ta denotes the time 
with beam on target reduced by the data acquisitions dead 
time (8 and 17 % for the lowest and highest explored en- 
ergies, respectively) and is known with an error of ±5%. 
edei stands for the detection efficiency which takes into ac- 
count the DSSSD detection efficiency (96.2 ± 0.1%) ||9], 
the ^^^Si*^^ charge state fraction which corresponds at both 
energies to 35% and finally the acceptance and the trans- 
port efficiency of DRAGON. This third parameter can be 
extracted using the GEANT simulations with the best agree- 
ment with data and corresponds to 32% (resp. 36%) for 
Ec.m. - 6.6 MeV (resp. 7.2 MeV). The error on the accep- 
tance and the transport through DRAGON is of the order 
of 15%. 
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Fig. 7. RC cross-sections (blue filled circles) along with fu- 
sion cross-sections (red crosses) LIU . Error bars for the fusion 
cross section correspond to the size of the red crosses. RC cross- 
sections above Vg are from a previous study (5). 



The obtained cross-sections are given in Fig. |7] along 
with data from our previous study [5] above Yb {Ec.m - 
8.5, 8.8 and 9.0 MeV) and fusion cross-section at these en- 
ergies extracted from 1 1 1 1. For the two recently explored 
energies RC cross-sections are lower than 1 ^b. Further- 
more excitation function slopes show that RC cross-section 
tends to decrease faster than the fusion cross-section below 
Vb- Indeed at Ec.m.- 9.0 MeV RC cross-section represents 



10 X 10 ^ of the fusion cross-section, this ratio falls to 
2 X 10"^ for the two energies discussed in this paper. 



4 Conclusion 

We have performed a heavy ion radiative capture reaction 
between two light heavy ions, '^C and '^O, leading to ^^Si 
at two resonant energies below Vb- Obtained y-spectra dis- 
play a previously unobserved strong feeding of interme- 
diate states around 11 MeV at these energies. This new 
decay branch has to be compared to the feeding of the 
same region observed at higher energies, in order to ob- 
tain the evolution of the direct feeding of the different ^^Si 
known states from the entrance-channel. As the new de- 
cay branch can not be fully reproduced by statistical nor 
semi-statistical decay scenarii this may be the signature of 
structural effects. With the upcoming new scintillator gen- 
eration, such as the PARIS project lfT2l . a sufficient res- 
olution will be achieved in order to disentangle the feed- 
ing of intermediate states around 1 1 MeV contained in the 
large bump. New Monte-Carlo simulations with other con- 
ditions such as calculated branching ratios with a cluster 
model for the entrance-channel have to be tested to see if a 
better agreement with the data can be achieved. 

GEANT simulations have been used to extract the DRA- 
GON spectrometer acceptance which is crucial to obtain a 
value for the radiative capture cross-sections. If new decay 
scenarii can reach better agreement with data, error bars 
on capture cross-section will be reduced. This will lead to 
a better understanding of the fusion reaction below Vb and, 
more particularly, for astrophysical energies where new 
resonances have been observed mainly in the neighbour- 
ing '^C-i-'^C system in the vicinity of the Gamow window 

ma. 
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